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ABSTRACT 


This  report  considers  the  problem  of  steady,  axial  tr  Nation  of  rigid 
and  fluid  spheres  in  stationary  and  moving  viscous,  incomMres  Jble  fluids 
bounded  by  an  infinitely  long  cylinder.  The  investigation  is  based  on  Stokes* 
approximation  for  the  hydrodynamic  equations  for  slow  flow;  thus  inertia  terms 
can  be  neglected,  and  the  stream  function  satisfies  a  fourth-order  differential 
equation  similar  in  form  to  the  bi harmonic  one. 

An  exact  solution  for  the  motion  of  rigid  spheres  m  s  ill  and  moving 
liquids  within  a  cylindrical  container  has  been  obtained  in  terns  of  an  infinite 
Rv-t  of  linear  algebraic  equations  for  the  coefficients  in  the  Stokes  stream  func¬ 
tion.  It  is  shown  that  the  drag  of  a  sphere  in  motion  within  a  moving  liquid  is 
composed  of  two  parts:  namely,  the  drag  due  to  the  motion  ot  the  sphere  ia  a 
still  Hauid  inside  the  cylindrical  tube,  and  the  drag  due  to  the  motion  of  the 
liquid  within  the  cylindrical  tube  past  a  stationary  sphere.  The  drag  experienced 
ij-  the  rigid  spheres  has  been  determined  for  the  two  special  cases  over  a  range 
of  ratios  of  sphere-to-cylindvr  diameter.  It  is  also  shown  that  the  first  two  equa- 
tte'i.r  of  the  infinite  set  closely  approximate  the  motion  a&tl  of  rigid  spheres 
over  a  larpe  range  of  diameter  ratios. 

For  fluid  spheres  (i.e.,  spheres  which  have  diffvic^i,  physical  properties 
than  the  external  medium  and  are  characterised  by  internal  motion)  an  approxi¬ 
mate  solution  (similar  to  the  one  for  the  rigid  case)  is  obtained.  The  drag  ex¬ 
perienced  by  the  fluid  spheres  has  been  computed  for  the  two  special  cases 
mentioned  above  for  a  range  of  diameter  ratios. 

Experimental  results  for  the  rigid  and  fluid  case  confirm  the  theory.  In 
general,  the  results  show  that  the  wall  effect  for  fluid  spheres  is  less  than  for 
corresponding  rigid  spheres.  Streamlines  and  velocity .drjrfr'Sr .tions  for  several 
cases  where  the  diameter  ratio  is  0.5  are  compared  with  those  in  an  infinite 
medium. 


INTRODUCTION 

The  effect  of  the  proximity  of  the  container  walls  on  the  drag  of  moving  bodies  is  of 
interest  in  many  fields  of  physics  and  engineering.  Examples  are  the  aae  of  the  failing  ball 
viscosimeter,  and  the  rise  of  air  bubbles  in  tubes.  The  purpose  of  the  present  investigation 
is  to  determine  the  drag  of  a  spherical  body  rising  or  falling  in  steady,  slow  motion  (i.e.,  the 
inertia  terms  in  the  equation  of  motion  can  be  neglected)  in  a  liquid  (stationary  or  moving) 
inside  an  infinitely  long  cylinder.  Both  rigid  and  fluid  spheres  a*r  bubbles)  are 

considered. 


? 


Previous  investigations  on  the  effect  of  the  proximity  of  the  container  walls  on  the 
drag  of  moving  bodies  have  dealt  with  rigid  spheres. •  Ladenburg1  and  Faxdn2  studied  the 
drag  of  spheres  moving  in  a  still  liquid  contained  in  an  infinitely  long  cylinder.**  Wakiya5 
and  Happel  and  Byrne6  determined  the  drag  of  spheres  in  Poiseuille  flow  (parabolic  velocity 
distribution)  in  a  cylindrical  tube.  Wakiya  considered  the  wall  effect  of  a  fixed  rigid  sphere 
in  Poiseuille  flow.  Happel  and  Byrne  also  included  the  case  where  the  rigid  sphere  is  moving 
inside  the  cylinder.  All  used  the  method  of  reflection  to  obtain  their  solution.  Starting  with 
the  known  solution  for  the  rigid  sphere  in  an  infinite  medium,  a  “reflection”  flow  is  superposed 
such  that  the  boundary  conditions  on  the  cylinder  are  satisfied  exactly.  The  boundary  condi¬ 
tions  on  the  sphere  are  only  approximately  satisfied.  The  drag  of  the  spheres  is  obtained  from 
Stokes’  law  using  the  velocity  of  the  sphere  increased  by  the  average  “reflection”  velocity 
on  the  sphere.  In  all  instances  approximate  expressions  for  the  drag  of  the  rigid  spheres  were 
given.  Cunningham,7  Williams,*  and  Lee9  considered  the  motion  of  a  rigid  sphere  at  the  in-  ’ 
stant  it  passes  the  center  of  a  spherical  container.  Bond10  suggested  an  approximate  expres¬ 
sion  for  the  wall  effect  of  fluid  spheres  in  cylinders,  based  on  LadenburgVs  results. 

In  the  present  study,  two  problems  are  considered.  The  first  problem  deals  with  the 
steady,  slow  motion  of  rigid  spheres  along  the  axis  of  an  infinitely  long  circular  cylinder. 

The  boundary  conditions  in  this  case  are:  uniform  velocity 'on  the  surface  of  the  sphere,  zero 
velocity  on  the  cylinder,  zero  velocity  or  parabolic  velocity  distribution  at  infinity.  The  solu¬ 
tion  of  the  problem  is  effected  by  means  of  the  Stokes  stream  function  for  axisymmetric  flow. 

Th©  motion  of  the  spheres  in  a  cylinder  is  solved  in  terms  of  a  system  of  linear  algebraic  equa¬ 
tions  for  the  constants  in  a  series  expansion  for  the  stream  function.  Numerical  results  are 
obtained  for  the  drag  of  the  rigid  spheres  for  diameter  ratios  up  to  0.8.  Values  of  wall  correc¬ 
tion  factors  have  been  computed  for  two  special  cases:  namely,  (1)  when  the  sphere  is  in 
motion  in  a  stationary  liquid,  and  (2)  when  the  liquid  is  in  motion  within  the  cylindrical  tube 
past  a  stationary  sphere.  The  drag  of  moving  spheres  within  a  moving  liquid  can  be  obtained  by 
appropriate  combination  of  these  two  correction  factors.  It  is  also  shown  that  the  first  two 
equations  of  the  infinite  set  closely  approx  mate  the  drag  over  a  large  range  of  diameter 
ratios. 

The  second  problem  deals  with  the  mot. on  of  fluid  spheres  (i.e.,  spheres  which  con¬ 
sist  of  a  fluid  that  has  different  physical  properties  than  the  external  medium,  and  are 


•During  publication  of  thi*  a  papar  daaltar  with  th*  wall  affact  of  ipturoid*  baa  coa*  to  tha  attaatioo 

of  tha  authors  (K^iVwofiCO  30), 

••Tha  comt^atHlinf  two-dime  o*?<5naI  problam  (a  circular  cyltedar  moo  tag  par  alia  I  with  tha  two  cootaiaar  walla) 
has  boon  traatad  by  Fa***2  sm.*  TaHut.i  * 

^Fafaroaccs  ar«  li*i*4  paga  65, 
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characterized  by  internal  motion)  along  the  axis  of  a  circular  cylinder.*  In  contrast  to  the 
rigid  case,  the  shape  of  fluid  bodies  ~».f\not  be  specified  beforehand.  The  shape  is  a  conse¬ 
quence  of  the  motion  such  that  the  boundary  conditions  are  satisfied  at  the  body  interface. 

For  the  slow  motion  of  a  fluid  body  in  an  infinite  medium,  it  was  shown  analytically  by 
Hadamard12  and  Rybczynski 13  that  the  sphere  is  a  possible  shape  for  which  the  boundary  con¬ 
ditions  (continuity  of  velocity  and  stress,  no  diffusion)  are  satisfied.  Experimental  evidence 
by  Spells14  indicates  that  the  spherical  shape  is  actually  attained  by  fluid  bodies  in  slow 
motion.  In  the  present  analysis,  a  spherical  shape  is  assumed  for  fluid  bodies  moving  within 
a  cylindrical  tube,  and  an  approximate  solution  (similar  to  the  approximate  solution  for  the 
rigid  case)  is  obtained.  As  before,  values  of  wall  correction  factors  have  been  computed  for 
the  two  special  cases.  The  validity  of  the  approximation  is  confirmed  by  experiments  with 
drops  moving  in  a  still  liquid.  The  experimental  investigation  shows  that  the  shape  of  a  fluid 
body  in  slow  motion  inside  a  cylinder  is  not  exactly  spherical.  The  drag  of  such  nonspherical 
bodies  can  be  evaluated  from  the  theoretical  solution  if  an  equivalent  radius  based  on  its 
volume  is  used,  and  if  the  diameter  ratio  is  smaller  than  0.5.  In  general,  it  is  shown  that  the 
wall  effect  for  fluid  spheres  is  less  than  for  corresponding  rigid  spheres.  Stream  functions 
have  been  evaluated  for  rigid  and  fluid  spheres  in  a  stationary  liquid  for  a  diameter  ratio  of 
O.b  and  are  presented  in  the  form  of  streamlines  about  the  spheres. 

The  investigations  described  in  this  report  were  carried  out  at  the  David  Taylor  Model 

Basin  under  NS  715-102  in  connection  with  a  program  investigating  gas-bubble  dynamics.  The 

|  ^ 

study  dealing  with  the  motion  of  spheres  in  a  still  liquid  was  first  presented  in  thesis  form 
and  subsequently  in  abbreviated  form.16  The  results  pertaining  to  the  motion  of  spheres 

a  m 

within  a  moving  liquid  have  been  presented  in  abbreviated  form. 

METHOD  OF  SOLUTION 

The  motions  considered  here  are  assuaged  to  be  sufficiently  slow  so  that  the  Stokes 
equations  of  motion  are  an  accurate  description  of  the  flow.  For  the  steady,  axisymmetric 
motion  of  an  incompressible  fluid,  the  Stokes  stream  function  exists  and  satisfies  a  fourth' 
order  differential  equation.  Corresponding  to  the  two  boundary  shapes  (cylinder  and  sphere), 
solutions  for  the  stream  function  are  obtained  in  cylindrical  and  spherical  coordinates.  In 
cylindrical  coordinates,  the  dUVential  equation  for  the  stream  function  is: 


•It  is  known  that  drops  sod  bubbles  below  cert* in  critics!  sixes  depart  frow  being  " fluid”  and  become  “rigid” 
in  thwir  behavior.  (Sea,  e.g..  the  am  bar  s’  investigation  on  the  motion  of  go*  bubbles.  Reference  II.)  It  is  not  the 
pwpoos  of  tbs  present  Investigation  to  dool  with  this  transition  phenomenon  nor  to  dotomino  the  conditions  ot 
which  transition  to  “rigidity”  totes  place.  It  is  intended  to  give  the  well  effect  for  those  droplets  or  bothies  that 
behave  as  fluid  bodies.  For  those  behaving  as  rigid  spheres,  the  well  effect  for  the  rigid  cooe  applies. 
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It  will  bo  shown  that  a  solution  of  this  equation  (for  the  case  of  symmetry  about  the  horizontal 
y-i  plane)  is: 


¥  (ar,  p)  •  ^  [p  K  j  (op)  f  (a)  +  p2  Kq  (op)  F  (o)  +  p  /j  (op)  (o)  4-  p2  /Q  (op)  G  (o)3  cos  axda  [2] 


Li  spherical  coordinates,  the  stream  function  satisfies  the  differential  equation 
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a  solution  of  which  is  shown  to  be 
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where  x,  p;  r,  a  are  coordinates, 


FI  ri  fi  a 

V  W  i 

f(c).F(a)>g(a),Cr(a) 


is  Gegenbauer  polynomial,*  order  n,  degree 
are  constants, 

are  modified  Bessel  functions,  and 
are  arbitrary  functions. 


For  the  satisfaction  of  the  boundary  conditions  on  the  cylinder  walls,  the  cylindrical  coordi¬ 
nate  solution  for  the  stream  function  is  used.  The  expression  thus  obtained  represents  the 
flow  inside  a  circular  cylinder,  not  as  yet  fully  specified  but  satisfying  the  boundary  condi¬ 
tions  on  the  cylinder.  This  expression  is  then  transformed  into  spherical  coordinates.  By 
comparing  termwise  the  constants  in  the  above  expression  with  the  constants  in  the  stream 
function  expansion  obtained  directly  in  spherical  coordinates,  a  relationship  between  the  con¬ 
stants  is  obtained.  The  boundary  conditions  on  the  sphere  yield  a  relationship  between  the 
constants  h  the  spherical  coordinate  solution.  Substituting  the  previous  relationships  into 
those  obtained  from  the  boundary  conditions  on  ihe  sphere,  a  set  of  linear  algebraic  equations 
'for  evaluating  the  constants  is  obtained. 


•TIm  first  tn  of  thm  ptlywlrto  Wm  •vaioata*  la  Ap»aadia  A. 
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THE  STREAM  FUNCTION  IN  SPHERICAL  COORDINATES 


The  Stokes  stream  function  in  spherical  coordinates  (Figure  1)  satisfies  the  fourth- 
order  differential  equation  (Reference  18,  p.  393): 

9  9 2 
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A  solution  of  the  differential  equation  is  sought  in  the  form 


(33 


V-l'i+Tj  (53 

such  that  L  Vj  *  0  and  L  Vj  *  1fl  (i.e.,  L2  V2  *  0) 

("a2  i-<2  a2 

where  L  *  I -  4- -  - 

[dr2  r2  dt2 

t  *  cos  0. 

The  solution  is  to  be  single-valued  and  continuous  in  the  given  Dow  field.  Assume  a  solution 
for  LVj  ■  0  exists  in  product  form;  i.e., 

? jCr^-JNr)  rU)  [6] 

The  second-order  partial  differential  equatioo  can  then  be  separated  into  two  ordinary  differ¬ 
ential  equations 


r2R"-a2R.  0  [73 

(<2  -  1)  f  a2  F  «  0  (83 

where  the  primes  denote  differentiation  with 
respect  to  the  argument. 

The  solution  of  Equation  [73  is: 

ft  (r)  -  Cj  r*  4-  c2  r~*+l  [9] 

1  >/l  +  4tf2 

where  *  *  -  ♦  - 2 - 


Figure  1  -  Coordinate  Systems 


The  solution  of  Equation  [8]  is:* 


Tti)  *C~*<J)  4  c3ra(t) 


I 


where  C~-  -  Gegenbauer  function  of  order  o,  degree  -*£,  and 
a2  *  «(o-  1), 

Comparing  the  expressions  for  m  and  a,  we  obtain  m  .  Hence 


Vx  if,  $)  -  (ct  r®4  c2  r-"+l)  (<?-*  4-  c3 


[ID 


The  solution  of  the  equation  L  ¥2  ~  ¥  j  is 


¥, *  V.  4  V 
2  1  P 
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Assume  the  particular  solution  ¥p  in  the  form  S(r)  T(t ),  hence 


and 


Therefore 
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Summing  the  above  solution  and  using  i  nteger  values  of  a  ,  the  expression  for  the  stream  func¬ 
tion  for  axisvmmetrie,  slow  flow  in  spherical  coordinates  becomes; 
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npressio;  $.  for  the  radial  and  tangential  velocities  of  the  flow  field  are  obtained: 
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[22] 

[50] 

[51] 

[52] 

Since  the 

im 


c,c3  *  gx(a),  c,  c4  -  fx 4a),  c%  c%  ~  Gx(a)t  cxc6  *  Fj(a),  c2  e3  -  ^  (a), 
c2  c4  "  r2cs  *  Gjfa)*  »nd  r2 cf>  •*  F 2(/)'  **l***sion  for  the  stream  function  in  cy¬ 

lindrical  coordi nates  becomes: 

T(x,  p)  -  I  [p  K ,  5  j  la)  4  p2  ^  F. '  )  -  p/j  <7j  (a)  4  p2  /0  («)]  cob  ax  da 

•* o 

+  J  [p ATj  /2 (a)  4  p2KQF2(a)  4  p/j  ?2  (o>  4  p2/0t?2(o)3  sin  axrfa  134} 


+  /fej  p2  4  Jt2  p4  4  kyp2  X  4  k4  p4X  4  4r5*2  4  k6p2X2  +kyX2  4  4g  * 

The  terms  containing  /(a)  and  p(a)  satisfy  the  second-order  differential  equation  (L  ?  -  0); 
the  F(o)  and  G(o)  terms  satisfy  the  fourth-order  differential  equation  (t2f  -  0).  The  arbitrary 
functions  /(a),  F(a),  y(a),  and  G(a)  are  evaluated  from  the  boundary  conditions  of  the  problem. 
From  the  relations 

I  — 

*  p  dp 


(35] 


1 

P  ds 


the  velocities  in  the  axial  and  p-direction  are  obtained: 


4 

-  J  [oK0fx(a)  •  {2K0  -  apKx)  Fj(o)  -  o10gtia)  —  (2  /Q  4  op/j)  Gj(«)l  coe  ax  da 
•*o 

4 1  Utf0/2(o)  - (2«f0-  eplfj)  F2(«)  -  al0g2(a)  - (2/0  4  op/j)  (?2(o)]  sin  a*  do 


-  2*j  -  4 *2p2  - 4, ■?  *  * Ff,  ~"z  ~  x2 

-  j  (Ay^cl  +  p^Fjt  ^  '  v?t*  >  ♦  pl0Gxi*))  asm 
•'o 

r„ 


/2<«)  ♦plT0Fj(^^f1f2(^  rp/0(72(o)U  cos  axis 


The  expression  for  the  pressure  field  is: 


[l0Gl(a)  4  h'Q  Fx  («)]  asinaxda 


2  ft  I  [IqG2Ui)  4  h'Q  F2  (o)l  ocos  at  da  4  p  (-16*2  x-%k4  x2  +ik*p2 -ik6z  46*7lnp) 
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THE  MOTION  OF  SPHERES  IN  A  CYLINDRICAL  TUBE 

A.  MOTION  IN  A  STATIONARY  LIQUID 
1.  Rigid  Spheres 

o.  Exoct  Solution.  The  first  problem  considered  is  the  motion  of  a  sphere  within  a 
still  liquid  of  finite  extent,  as  occurs  when  a  sphere  rises  or  falls  under  the  influence  of 
gravity.  In  this  section  the  axial  motion  of  a  rigid  sphere  inside  an  infinitely  long  cylinder 
will  be  dealt  with.9 

For  convenience,  the  coordinate  origin  is  taken  at  the  center  of  the  sphere  and  the 
cylinder  is  assumed  to  be  moving  at  constant  velocity  V  in  the  negative  e-direction  (Figure  2). 
The  boundary  conditions  are: 

at  the  cylinder  walls  (p  *  A):  v  «  -£/;  v  «  0  or  ¥  «—  b2 

x  P  2 

at  infinity  (*  «  ±m) •  e  «  -U\  t>  -  0  [381 

*  P 

on  the  surface  of  the  sphere  (r  *  H)i  vf  *  0;  e^  «  0 

The  case  of  a  sphere  moving  at  constant  velocity  U  is  obtained  bom  the  above  by  superimpos¬ 
ing  a  uniform  flow  with  velocity  U  in  the  positive  e-direction. 

Since  the  stream  function  is  symmetrical  about  the  y- »  plane  because  of  symmetry  in 
the  boundary  conditions^  the  stream  function  in  spherical  coordinates  is  given  by 
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Sphere  Moving  at  Velocity  V  In  Still  Liquid 
In  aid*  Fixed  Cylindrical  Tuba. 


Equivalent  System:  Sphere  Fixed;  Tube  Moving 
at  Velocity  U.  Uniform  Velocity  V  at  Infinity. 


Figure  2  —  Definition  Sketch 


The  boundary  conditions  on  the  cylinder  walls  can  be  satisfied  if  we  take  the  stream  function 
in  cylindrical  coordinates  as 


f  ix,  p) 


[p  K  t  (ap)  fx  (a)  ♦  p2  K0  («p)  Fx  (o) + p  t  y(ap)  g ,  (a)  +  p7  l0iap)Gl  (a)]  cos  as  da  ♦ 

1401* 


Utilising  the  boundary  conditions  on  the  cylinder  walls,  (Equation  lS8l),  we  obtain  the  follow¬ 
ing  relations: 

(ram  T  «  -»J  :  »*,(«*)/,(«)  +  »*  *„(«»)/■,(•)  +  »/,(«*)*,  to  ♦  ^UWGjto -0 
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If 
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10  l 


W'e  can,  therefore,  solve  for  gx  (a)  and  G  { (a)  in  terms  of  fx  (a)  end  Fx  (a) 
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where  the  following  relations  for  the  Bessel  functions  have  been  used: 
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Thus  having  satisfied  the  boundary  conditions  on  the  cylinder  walls,  the  stream  function 
becomes 


+  p2  l0(ap)[S2fl(a)  +  S4Fx  (a)]}  cos  ax  da  4  ~ 


In  order  to  transform  the  above  expression  into  a  form  amenable  to  the  satisfaction  of  the 
boundary  conditions  on  the  sphere,  let 
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The  integrals  prf  KAap)a*  cos  ox  da  and  p  f  K .  {a p)  a*  cos  ax  da  which  then  resul  t  are 

evaluated  in  Appendix  B.  It  is  shown  there  that  only  odd  powers  of  a  for  ^  and  even  powers 
of  a  for  bm  are  of  interest  here;  otherwise  discontinuities  result  at  $  «  0,  #  Thus  wo  obtain 
from  Equation  [40] 
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Ip2  K  Q(  *p)  ($0  +  o2  -4-  ft4o4  4  .  .  .  )  cos  as  da 
o 

1  [plxUp)gx 


1473 


(o)  4  p2  l0(ap)  G  x  (o)l  cos  ax  da  4 
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The  integrals  containing  /Q  and  can  be  expanded  into  a  Taylor  series  about  the  origin.  The 
convergence  of  these  integrals  is  assured  by  the  form  of  the  functions  0j(a\  and  C x  (a).4 
A  function  ^  (x,  p)  can  be  expanded  by  Taylor's  theorem  into 
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is  thus  transformed  into 
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From  the  solution  in  sohericai  coordinates  (equation  (*•],  we  obtain  (making  use  of  the  relations 
p  ft  *  sin  $  and  x/r  «  cos  0): 
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A  number  of  the  S  *  and  S*  integrals  have  been  evaluated  in  Appendix  C  snd  are  tabulated 
in  Table  3.  Thus  we  obtain 
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From  the  boundary  conditions  on  the  sphere  (Equation  [881)  and  the  orthogonality  relations 
of  Gegenoauer  and  Legendre  polynomials,*  we  obtain 
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Solving  for  A  and  C%  ,  we  obtain 
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Substituting  the  expressions  for  B>,  C>t  D*  as  given  in  Equations  [58]  into  [55],  an 
infinite  set  of  linear  algebraic  equations  is  obtained  for  determining  the  constants  at,  «3, 
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The  constants  a  and  b  can  be  evaluated  numerically  to  any  degree  of  accuracy  by  increasing 

■w 

the  number  of  equations  utilized. 

The  solution  for  the  stream  function  in  spherical  coordinates  for  flow  about  a  sphere 
m  an  infinitely  long  cylinder  is  thus  given  by  Equation  [39]  with  the  constants  4#,  Bm>  Cm, 
and  Dm  given  in  Equations  [53].  Using  ten  equations  of  the  algebraic  set  (Equations  [56]),  a 
sufficient  number  of  the  constants  am  and  5#  was  determined  for  a  diameter  ratio  of  0.5. 
Streamlines  were  then  evaluated  from  Equation  [39].  They  are  shown  in  Figures  3  and  4. 
Figure  3  shows  the  relative  motion  about  the  sphere;  i.e.,  the  streamlines  as  they  appear  to 
an  observer  moving  with  the  sphere.  Figure  4  shows  the  absolute  motion  about  the  sphere; 
i.e.,  the  streamlines  as  they  appear  to  an  observer  fixed  in  space.  Velocity  distributions 
at  the  plane  of  symmetry  and  streamlines  for  the  motion  of  a  corresponding  sphere  in  an  infi¬ 
nite  medium  have  been  included  in  both  figures.* 


•TW  «sprMilM  tor  tko  *ro—  fuactioo  tor  tho 


otrifid 


?(r, *)»—  ola2*  (hr*+L  UR3  -  -  1  URr) 
2  V  2  r  2  / 

’P(r.d)*-  UR3  i-  — t/gr) 

2  \  2  r  2  / 


la  aa  Ufloitt  oMdioa  io 


rrltUrt  motlom 


20 


*** 


Sphere  in  Circular  Cylinder,  Diameter  Ratio:  0.5  Sphere  in  Infinite  Medium 

Figure  3  -  Streamlines  and  Velocity  Distribution  for  Rigid  Spheres 
(Exact  Solution);  Relative  Motion 


The  drag  experienced  by  moving  spheres  can  be  evaluated  by  integrating  the  forces 
meting  over  the  surface  of  the  sphere.  Thus 


Tangential  Drag 


mhere  prg  is  the  tangential  stress, 
»  .  is  the  normal  stress,  and 

fr 

dS  is  the  surface  element. 


"The  drag  of  the  spheve  can  be  obtained  from  the  stream  function  and  is  shown  in  ;  **ndix  0 
go  be  equal  to 

Drag  -  [58] 


i.e.  *  the  drag  is  proportional  to  the  coefficient  02  of  the  stream  function  expansion.  Defining 
s  wall  'ejection  factor  (Kx)  as  the  ratio  of  the  drag  of  the  sphere  in  the  bounded  medium  to 
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Figure  4  -  Streamline*  and  Velocity  Distribution  for  Rigid  Spheres 

(Exact  Solution);  Absolute  Motion 


that  in  an  infinite  medium,  the  drag  can  be  expressed  as 
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Also,  combining  Equations  (5S]  and 
factor  in  terms  of  the  coefficient  bQ: 
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f  we  obtain  an  expression  for  the  wall  correction 
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Wall  correction  factors  for  rigid  spheres  moving  in  a  still  liquid  inside  an  infinitely  long  cyl¬ 
inder  have  been  determined  by  numerically  solving  the  algebraic  system  (Equation  l5ij)  for 
the  coefficient  bQ  over  a  large  range  of  diameter  ratios*  The  number  of  equations  of  the  alge¬ 
braic  system  used  was  increased  (at  most  up  N  eijrt)  vrtfl  only  very  small  changes  in  the 
value  of  t0  were  obtained.  The  oomputed  wall  ccnecffon  Actors  are  shown  is  Figure  S  and 
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TABLE  1 


Wall  Correction  Factors  (Kt)  for  Rigid  Spheres  Moving  in  «  Still  Liquid 

inside  s  Cylindrical  Tube 


A 

Exact  Theory 

Approximate  Theory 
(Equation  162]) 

Percent  Difference 

0.0 

1.000 

1.000 

0.00 

0.1 

1.263 

1.263 

0.00 

0.2 

1.680 

1.680 

0.00 

0.3 

2.371 

2.370 

-  0.04 

0.4 

tm 

3.582 

-  0.39 

0.5 

5.970 

5.871 

-  1.66 

0.6 

11.135 

10.591 

-  4.89 

0.7 

24.955 

21.40C 

-14.22 

0.8 

73.555 

_ _ _ 

48.985 

-33.40 

Table  1.  Previous  theoretical  results  by  Laden  burg  and  Faxdn  are  included  in  Figure  5.  Ex¬ 
perimental  data  at  very  low  Reynolds  numbers,  taken  from  references  9,  22,  23,  and  24.  are 
also  shown  for  comparison.  Only  experimental  data  at  Reynolds  number  (based  on  diameter) 
less  than  two  were  included.  The  actual  extent  of  the  Reynolds  number  range  over  which  the 
theoretical  solutions  for  the  drag  of  the  spheres  are  valid  could  not  be  determined  from  avail¬ 
able  experimental  data.  It  is  estimated,  however,  that  the  theoretical  solutions  are  reliable 
for  Reynolds  numbers  up  to  two. 


b.  Approximate  Solution.  A  very  good  approximation  for  the  drag  of  the  rigid  spheres 
can  be  obtained  by  retaining  only  the  first  two  equations  of  the  infinite  set  (Equation  [56]) 
and  considering  only  the  constants  b  and  o  j .  For  the  constant  60  we  then  obtain 
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Using  the  definition  of  (Equation  [60]),  the  wall  correction  factor  lor  rigid  spheres  in  s 
still  liquid  within  an  infinitely  long  cylinder  is  then  given  lathe  coaveaient  form: 
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Wall  correction  factors  given  by  Equation  [62]  have  been  evaluated  and  are  included  in  Fig¬ 
ure  5.  In  Table  1,  wall  correction  factors  for  the  exact  and  approximate  solution  are  listed  for 
comparison.  It  is  seen  that  the  wall  correction  factors  for  rigid  spheres  is  well  approximated 
by  Equation  [621  up  to  diameter  ratios  of  0.6. 

In  this  approximation  the  stream  function  becomes 
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In  the  solution  given  in  Equation  [691,  all  boundary  conditions  on  the  cylinder  are  satisfied. 
The  velocity  components  arising  from  the  C~*  -terms  satisfy  the  boundary  conditions  on  the 
sphere.  Those  from  the  higher  terms  (C  ...  ,  etc.)  do  not  satisfy  the  boundary  conditions 
on  the  sphere. 

Streamlines  were  evaluated  from  Equation  [69]  for  a  diameter  ratio  of  0.5,  and  are  shown 
in  Figure  6  for  relative  motion  about  the  sphere;  i.e.,  as  they  appear  to  an  observer  moving 
with  the  sphere. 
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Sphere  in  Circular  Cylinder, Dimeter  Ratio;  0.5  CMererieee  of  Velocity  Dietrlbatioe 

Figure  9  -  Streamlines  and  Velocity  Distribution  fer  Rigid  Spheres 
(Approximate  Solution);  Relative  Motion 


2.  Fluid  Spheres 

In  this  section  the  axial  motion  of  a  fluid  body  inside  an  infinitely  long  cylinder  wilt  be 
considered.  For  the  motion  in  aa  infinite  medium,  it  has  been  shown  that  the  sphere  represents 
a  possible  snrface  at  which  all  boundary  conditions  are  satisfied.1*’15  It  will  be  shown  hers 
that  the  assumption  of  a  spherical  shape  doee  not  lend  to  an  enact  solution  for  the  motion  of 
a  fluid  body  inside  n  cylinder.  However,  nil  velocity  and  stress  compos  sals  srisiag  from  the 
first  group  of  terms  of  the  stream  fuactioa  (the  Cf*  terms)  do  satisfy  the  bouadary  coaditioa* 
on  a  sphere.  Hence,  for  the  case  when  the  Cf*  terms  see  most  important  in  describing  the 
motion,  the  spherical  shape  will  be  a  very  good  approximation.  U  approximate  solntion  similar 
to  the  one  for  the  rigid  case  will  thus  be  obtained*  Tbs  validity  of  this  eolation  is  confirmed 

by  experiments. 
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*  P 

at  the  surface  of  the  fluid  body:  (6 1 1 

Equality  of  tangential  velocities, 

Normal  velocities  vanish  (i.e.,  no  diffusion), 

y.  1 

Equality  of  tangential  and  normal  forces. 

•  K 

For  a  spherical  surface  (r  -  R),  the  boundary  conditions  of  Equatin’.  ^  He  frrm: 
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where  uf  is  the  radial  velocity  inside  sphere, 

u^is  the  tangential  velocity  inside  sphere,  and 
i  refers  to  quantities  inside  sphere. 

The  stream  function  in  the  exterior  of  the  sphere  is  «  * 
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From  the  above  expressions  for  the  stream  function  at  csjres* 
poncnts  (Equations  [18])  and  stress  components  (Equation  -j*  and 
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obtained  at  the  Center  o#  the  sphere. 
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In  spherical  coordinates  it  is  given  by 

¥(r,  (9)  «  C  (cos  (?)  L2  f2  +  «2  ~  +  C2  f4  +  02  rj  ♦  C?4“*  (cos  0)  ^4  r4  ♦  C4  r*J  ♦ .  , 
From  the  boundary  conditions  on  the  sphere  (Equation  [681,  *  *  2) 
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where  a 


external  viscosity 
internal  viscosity 


Substituting  expressions  for  4a,  Ca,  Br  and  D%  (Equation  [«3])  into  Eqnotion  [70],  two  alge¬ 
braic  equations  for  ax  and  &0  arc  obtained 
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From  the  above  equations  &0,  for  example,  becomes 
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The  constants  bQ  and  av  as  determined  from  Equation  [71],  when  substituted  in  Equation  [63] 
will  yield  the  exterior  stream  function  for  the  fluid  sphere.  In  the  interior  of  the  sphere  the 
stream  function  can  be  given  as  (utilising  velocity  continuity): 


<t  (r,  0  -  C^‘  (cos  9)  Ea  f*  ^  1  -  -M  +  C'4~*  (cos  9)  M4  f4  +  CA  f5)  +  . . 


ITS] 


Streamlines  for  the  interior  and  exterior  of  the  fluid  sphere  were  evaluated  from  Equations  [73] 
and  [63]  (with  coefficients  bQ  and  of  determined  from  Equation  [71])  for  the  case  of  an  infinite 
viscosity  ratio  and  a  diameter  ratio  of  0.5.  They  are  shown  in  Figures  8  and  ?.  Figure  8 
dhows  the  relative  motion;  Figure  9  gives  the  absolute  motion.  Velocity  distributions  at  the 
plane  of  symmetry  and  streamlines  for  the  motion  of  a  corresponding  sphere  in  an  infinite 
medium  have  been  included  in  both  figures.* 

Using  Equations  [58]  and  [72],  we  obtain 
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Dr*f,  In  Infinite  Medium 


•The  street*  function  for  the  motion  of  fluid  ephetee  in  an  infinite  medium  is; 
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Sphere  in  Circular  Cylinder, Diameter  Ratio:  0.5  Sphere  in  Infinite  Medium 

Figure  8  -  Streamlines  and  Velocity  Distribution  for  Fluid  Spheres;  Relative  Motion 


where  Kl,  the  wall  correction  factor  for  fluid  spheres  moving  in  a  still  liquid  inside  an  infi¬ 
nitely  long  cylinder,  is 


a 


Wall  correction  factors  have  been  evaluated  from  Equation  [75 1  and  ire  give*  in  Figure  10. 
Curves  for  fluid  spheres  of  viscosity  ratios  0  (rigid  case),  1,  IS,  and  **  ere  shows. 


•For  a  *0.  tit  ttgHtatioft  for  K j  eteu  Is  that  af  Ipslw  (ill. 


SI 


Sphere  in  Circular  Cylinder, Diam*t»r  Ratio:  0.5  Infinite  Radium 

Figure  9  —  Streamlines  and  Velocity  Distribution  for  Fluid  bpheres;  Absolute  Motion 


b.  Experimental  Work.  To  check  the  validity  of  the  theoretical  solution  foe  the  drag 
of  fluid  spheres  and  to  determine  the  deviation  of  fluid  bodies  in  a  cylindrical  container  from 
the  spherical  shape,  a  limited  number  of  experiments  was  conducted.  The  experimental  study 
consisted  of  determining  the  rate  of  fall  and  the  site  of  water-glycerine  and  Dow-Coming  200 
silicone  drops  in  caator  oil.  The  DC  200  corresponded  to  viscosity  ratio  of  13,  whereas  the 
water-glycerine  mixture  (about  45  percent  glycerine)  gave  a  ratio  of  about  200.  Two  cylindri¬ 
cal  Lucite  tubes  were  used.  They  were  1.24  in.  and  2.74  In.  in  internal  diameter  and  36  in, 
and  31  in.  in  length,  respectively.  The  drop#  were  generated  by  means  ot  a  stopcock  burette, 
the  tip  of  which  provided  into  the  castor  oil.*  fly  regulating  the  etopcock,  drop*  of  different 
site  were  released.  Large  dropa  were  used  to  assure  that  they  behaved  as  “fluid  bodies,* 

The  fall  velocities  of  the  drops  were  determined  by  means  of  a  atop  walchj  the  site  and  shape 
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m 


df  the  drops  were  determined  from  high-speed  photographs.  Stationary  steel  spheres  of  various 
known  diameters,  located  at  the  center  of  the  tubes,  were  photographed  to  provide  scale  factors 
for  evaluating  sire  and  shape.  The  temperature  of  the  castor  oil  was  deterrd*l*d;  at  frequent 
intervals,  by  means  of  an  immersion  thermometer.  The  viscosity  of  the  liquids  was  measured 
with  an  Ostwald- Fen ske  viscometer,  and  the  density  of  the  \  inhb->  «  obtained  by  means  of 
a  Westphal  specific-gravity  balance. 

The  experimental  drag  of  the  drops  was  determined  from  the  equilibrium  condition  of 
the  forces  acting  on  them.  Thus 

Drag  *  Volume  yi  g  -  Volume  y  g 

\  y  ■  V—  .1-  — y  / 

Weight  Buoyancy 


The  experimentally  determined  drag  can  be  compared  with  the  theoretical  one  given  in  Equa¬ 
tion  [74].  To  obtain  the  experimental  wall  correction  factors,  the  experimental  drags  in  the 


cylindrical  tubes  were  divided  by  the  drag  of  the  fluid  spheres  in  an  infinite  medium.  For 


deformed  drops,  an  equivalent  radius 


was  used  in  the  computations. 


The 


experimental  correction  factors  thus  determined  are  shown  in  Figure  10.  The  drops  as  deter¬ 
mined  from  the  experiments  were  either  spheres  or  ellipsoids  of  revolution  with  major  axes  fn 
the  direction  of  motion.  The  eccentricity  of  the  drops  is  given  in  the  insert  of  Figure  10. 

The  maximum  Reynolds  number  (based  on  diameter)  of  the  water-glycerine  drops  was  0.72; 
th©  maximum  for  the  DC  200  drops  was  0.2.  It  is  seen  from  this  figure  that  the  theoretical 
solution  agrees  well  with  experimental  results  for  diameter  ratios  up  to  about  0.5.  The  devia¬ 


tion  from  the  theoretical  curves  occurs  at  diameter  ratio  of  0.53  for  the  water-glycerine  drops 


and  at  a  ratio  of  0.65  for  the  DC  200  drops.  Because  of  the  limited  range  of  the  experiments, 
no  definite  conclusion  can  be  drawn  regarding  this  difference.  It  is  conceivable,  however, 
that  the  deviation  from  the  approximate  theory  is  a  function  of  the  viscosity  ratio  a\  i.e.,  the 
smaller  this  ratio  the  larger  the  diameter  ratio  at  which  the  theory  deviates  from  the  experi¬ 
mental  results.  Although  the  largest  a  used  was  200,  it  is  not  believed  that  the  deformation 
of  bodies  wifh  a  -*  m  would  be  much  more  severe.  It  is  further  noticed  from  Figure  10  that  the 
wall  correction  factor  is  not  too  sensitive  to  deformation  from  the  spherical  shape.  A  differ¬ 
ence  in  the  axes  of  the  drops  of  up  to  about  15  percent  occurred  without  noticeably  affecting 
the  wall  correction  factor. 
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B.  MOTION  IN  A  MOVING  LIQUID 

The  preceding  sections  have  dealt  with  the  steady,  axial  translation  of  rigid  and  fluid 
spheres  in  a  stationary,  viscous,  incompressible  fluid  bounded  by  an  infinitely  long  cylinder. 
The  solutions  are  now  extended  to  include  the  case  in  which  the  fluid  contained  within  the 
cylinder  is  also  in  motion;  i.e.,  a  parabolic  velocity  distribution  (Poiseuille  flow)  exists  at 
infinity.  As  before,  two  problems  are  considered;  the  first  deals  with  the  motion  of  rigid 
spheres,  the  second  with  fluid  spheres.  The  method  of  solution  is  the  same  as  before;  the 
only  difference  arises  from  the  changed  boundary  conditions  at  infinity. 

An  exact  and  approximate  solution  is  given  for  the  drag  of  moving  rigid  spheres  in  a 
moving  liquid.  The  approximate  solution  is  shown  to  be  valid  for  diameter  ratios  up  to  about 
0.0.  An  approximate  solution  is  given  for  fluid  spheres  and  is  estimated  to  be  valid  for  diam¬ 
eter  ratios  up  to  about  0.5.  It  is  shown  that  the  drag  of  a  sphere  in  motion  within  a  moving 
liquid  is  composed  of  two  parts:  namely,  the  drag  due  to  the  motion  of  the  sphere  in  a  still 
liquid,  and  the  drag  due  to  the  motion  of  the  liquid  within  the  cylindrical  tube  past  a  station¬ 
ary  sphere.  Values  of  wall  correction  factors  have  thus  been  computed  for  the  two  special 
cases.  The  drag  of  moving  spheres  within  a  moving  liquid  can  be  obtained  by  appropriate 
combination  of  these  two  correction  factors. 


1.  Rigid  Spheres 


o.  Exoct  Solution.  The  problem  considered  here  is  that  of  a  rigid  sphere  moving  within 
a  fixed  circular  cylinder  containing  a  fluid  having  a  parabolic  velocity  distribution  at  infinity. 
For  convenience,  it  is  desired  to  take  the  coordinate  system  fixed  with  respect  to  the  sphere. 
Thus  an  equivalent  system,  as  indicated  in  Figure  11,  is  used.  The  coordinate  origin  is  again 


taken  at  the  center  of  the  sphere,  and  the  cylinder  is  assumed  to  be  moving  at  constant  veloc¬ 
ity  V  in  the  negative  sr-direction.  The  boundary  conditions  are: 


V  ,  V  , 
-0;  e  «0  or  *  *  _  ft2  -  -  62 

2  4 

2\ 


at  the  cylinder  walls  (p  *  b):  v 

*  i 

»t  infinity  (*  « t  *>):  +  K  (l  —  —\  ;  v  *  0 

\  w  p 

on  the  surface  of4he  sphere  (r  *  R):  vr  *0;  *  0 


[76] 


After  satisfying  the  boundary  conditions  on  the  cylinder  vails,  the  stream  function  becomes 
(similar  to  Equation  [45]) 


ipA  ,  (ap)  (a)  +  p 2  K0 (op)  F,  (a)  +  p  /,  (ap)  [S ,  /,  (a)  +  S, F,  (a) ) 


+  p1  /„  (ap)  IS,/,  (a)  +  S,  Fj  (a)]  |  oon  ax  da  +  ^ 


+  -  P 
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Sphere  Moving  with  Velocity  U  in  Moving  Liquid 
Inside  Fixed  Cylindrical  Tube 


Equivalent  System:  Sphere  Fixed;  Tube  Moving  with 
Velocity  U.  Velocity -U+  V  (1 ~p2/b2)  at  Infinity 


Figure  11  —  Definition  Sketch 


As  before,  the  stream  function  ¥  is  now  rewritten  as 
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Again,  comparison  with  the  solution  for  V  in  spherical  coordinates  (Equation  [52])  results  in 
the  following  relationships  between  the  constants  A  ,  B  ,  C  ,  D  ,  and  a  ,  b  : 
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etc. 


As  previously  (Equation  [55]) 

2n  +  1  1  2n  -  1  1 

A  **  —  - - B  - —  - - *  0  - 
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Substituting  the  expressions  for  Am ,  8  .  C  ,  0  as  given  in  Equation  [79]  into  [55],  the  infinite 
set  for  the  determination  of  the  constants  alt  <*3,  . . .  and  6Q,  62  ... ,  is  obtained 
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As  in  the  case  of  the  motion  of  a  sphere  in  a  still  liquid,  a  wall  correction  factor  could  be 
defined  based  on  the  drag  in  an  infinite  medium  [  i.e.,  6 » rp  R(U  -  V)].  However,  the  wall 
correction  factor  thus  obtained  would  be  a  function  of  the  diameter  ratio  X  and  the  velocity 
ratio  V/V.  It  is  more  convenient  to  define  a  correction  factor  which  is  a  function  of  X  only. 

The  coefficient  b0  as  obtained  from  Equations  [80]  can  be  written  (utilizing  properties 
of  determinants)  as: 

*o  -  \>i  -  *02 


where  6Q1  is  determined  from  System  1  and  bnj  from  System  2 
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-  (  -  1.07056  X3)  +  .  .  .  X2 
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As  given  previously,  the  drag  experienced  by  a  sphere  is 

Drag  *  -  4  *r‘  n  60 


wherr  the  drag  is  taken  positive  in  the  negative  ^direction  as  indicated  in  Figure  11.  We 
now  define  two  wall  correction  factors  Hx  and  k2  such  that 
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Walt  Correction  Fact**,  K 


TABLE  2 

Wall  Correction  Factors  (K2)  for  Fixed  Rigid  Spheres  in  Poiseuille  Flow 


- — 

X 

Exact  Theory 

Approximate  Theory 
(Equation  [16]) 

Percent  Difference 

0.0 

LOCO 

1.000 

0.00 

0.1 

1.255 

1.255 

0.00 

0.2 

1.635 

1.635 

0.00 

0.3 

2.23! 

2.231 

0.00 

0.4 

3.218 

3.218 

0.00 

0.5 

5.004 

4.973 

-  0.83 

0.6 

8.651 

8.377 

-  3.17 

0.7 

17.671 

15.666 

-11.25 

0.8 

47.301 

33.056 

-30.1 

hence 
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Figure  12  -  Wall  Correction  Factors  for 
Rigid  Spheres  within  a  Cylinder 
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Thus  the  drag  of  a  sphere  in  motion  within  a 
moving  liquid  is  composed  of  two  parts:  (1) 
the  drag  due  to  the  motion  of  the  sphere  with 
velocity  (J  in  a  still  liquid,  and  (2)  the  drag 
due  to  the  sphere  held  fixed  within  a  moving 
liquid  having  a  parabolic  velocity  distribution 
(maximum  velocity  K)  at  infinity. 

The  coefficients  bQl  have  been  deter¬ 
mined  in  v  previous  section  from  Equations 
[56].  The  coefficients  &02  were  determined 
numerically  from  System  2,  Equation  [81],  in 
the  same  manner  as  6Q1  over  a  range  of  diam¬ 
eter  ratios  X.  The  computed  wall  correction 
factors  A'j  (for  rigid  spheres  moving  in  a  still 
liquid  inside  an  infinitely  long  cylinder)  were 
shown  in  Table  1.  Table  2  gives  the  compute* 
wall  correction  factors  A'2  (for  fixed  rigid 
spheres  within  Poiseuille  flow).  Both  wall 
correction  factors  are  shown  in  Figure  12. 

The  stream  function  in  spherical  coor¬ 
dinates  for  flow  about  a  moving  rigid  sphere 
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In  Poiseuiile  flow  is  given  by  Equation  [30]  with  coefficients  A C%1  as  determined 
from  Equation  [79}  (with  coefficients  a*  and  evaluated  from  Equations  [80]  or  Bll), 


b.  Approximote  Solution.  Again,  an  approximate  solution  for  the  rigid  spheres  can  be 
obtained  by  retaining  only  the  coefficients  bQ  and  From  Equation  [80]  we  obtain 

{(J-V)  (1-  0 ,,75857  A5)  +  V[- 
£ _ _ \3 

R  2n  1-2.1050  A  +  2.0065  A2 3  -  17068  aVo.72603  A6 
The  drag  of  the  moving  rigid  spheres  in  a  moving  liquid  becomes 

W-V)  (1  -0.75857  As)  +  F^  A2  -0.55640 

Drag -6 n^R - — - _ - - -  * 

1-2.1050  A +2. 0665  A3 -1.7068  A5 +0.72603  A° 

losing  definitions  of  Kx  and  K2  (Equations  [60l  and  [82]),  we  obtain 

[62] 


[86] 


Wall  correction  factors  as  obtained  from  Equations  [62]  and  [86]  are  shown  in  Figure  12  and 
Tables  1  and  2.  Previous  theoretical  results  by  Happel  and  Byrne®  and  Wakiya5  (presented 
in  terms  of  K 2)  are  also  included  in  Figure  12. 

It  is  seen  from  the  tables  that  the  approximate  expressions  for  Ky  and  K2  agree  well 
with  the  exact  solution  up  to  diameter  ratios  of  0.6.  However,  for  certain  combinations  of 
€JK{  -VK2  (e.g.,  U/V  «  3/4),  the  range  of  good  agreement  will  decrease  to  diameter  ratios  of 

•bout  0.5. 

The  approximate  solution  for  the  stream  function  in  spherical  coordinates  for  flow  abou 
a  moving  rigid  sphere  in  Poiseuiile  flow  is  given  by  Equation  [68]  with  coefficients  Am,  Bm, 

0,0  as  follows: 
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mm  as  in  Equation  (63] 


2.  Fluid  Spheres 

For  fluid  bodies  moving  within  Poiseuille  flow,  the  same  type  of  Approximation  is 
utilised  as  for  the  motion  of  a  fluid  sphere  in  a  still  liquid.  As  before,  the  origin  of  the  coor¬ 
dinate  system  is  taken  at  the  center  of  the  fluid  body.  The  boundary  conditions  are: 

at  the  cylinder  walls  (p  -  b) :  v%  »  -  V ;  e  «  0 


at  infinity  (x  -  -«) : 


at  the  surface  of  the  fluid  body: 

Equality  of  tangential  velocities 

Normal  velocities  vanish  (i.e.,  no  diffusion) 

Equality  of  tangential  and  normal  forces. 


[88] 


Again  assuming  a  resultant  spherical  shape  for  the  fluid  body,  we  obtain,  as  before,  from  the 
boundary  conditions  on  the  sphere 
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Substituting  expressions  for  Av  Cy  and  02  (as  givon  in  Equation  [87])  into  Equat?D»>  f  1>- 
two  algebraic  equations  for  b0  and  ax  are  obtained 
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From  the  above  equations  &0  and  »  become 
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The  drag  of  the  moving  fluid  spheres  in  a  moving  liquid  then  becomes 
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the  wall  correction  factors  Kx  Mid  K  2  become 
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Figure  13  -  Wall  Correction  Factors  Figure  i  -  Wall  Correction  Factors  K2 

for  Spheres  within  a  Cylinder  for  Spheres  within  a  Cylinder 


[75] 


[94] 


thus 


Drag  «  6  wp.fi 


1  +  2/3  a 

1+  a 


uw,  -  vh2) 


[95] 


The  wall  correction  factors  Kt  and  k2  (as  obuined  from  Equations  [75]  and  [94])  are  given  as 
a  function  of  the  diameter  ratio  in  Figures  18  and  14  for  three  viscosity  ratios  (a  *  0  [rigid 
case],  1,  and  The  above  solutions  for  the  wall  correction  factors  for  the  fluid  spheres  are 
estimated  to  be  valid  for  diameter  ratios  up  to  0.5. 

The  expression  for  the  exterior  stream  function  for  fluid  spheres  in  Poiseuille  flow  is 
given  by  Equation  [87]  with  coefficients  6Q  and  at  determined  from  Equation  [89].  The  stream 
function  in  the  interior  region  is  given  by  Equation  [78]  (with  coefficients  as  given  in  Equation 
tt>7]  aad  determined  from  Equation  [St]). 
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TERMINAL  VELOCITIES  OF  SPHERES 

In  many  applications  it  is  of  interest  to  determine  the  terminal  velocity  (V)  of  rising  or 
falling  spheres.  Such  a  velocity  can  easily  be  determined  from  the  res  vita  of  the  present  in¬ 
vestigation.  A  body  rising  or  falling  under  the  influence  of  gravity  reaches  such  a  velocity 
(terminal  velocity)  when  all  forces  acting  on  it  are  in  equilibrium: 


or 


Drag  4-  Buoyant  Force  ♦  Weight  *  0 


6  WfiR  - -  WKy  -  VKJ  -  i  wR*  yg  +  -  nRZy.g  «  0 

1  +  <r  *  3  3  1 


m 


where  y  is  the  density  of  liquid  and  yi  is  the  density  of  sphere;  hence 


“  y4  ) 


l  +  <r 


[97] 


for  o  *  0,  V  *  0,  and  *  1;  i.e.,  rigid 
spheres  moving  in  a  still  infinite  medium, 
Equation  [97]  reduces  to  Stokes  law;  namely, 

2  R2  9 (y  -  y  ) 

V  - - - 

9  n 

The  value  of  the  appropriate  wall  correction 
factor  is  obtained  from  Figure  5  or  Table  1, 
Figure  10  and  Figure  14  or  Table  2.  The 
ratio  K2/Kx  appearing  in  Equation  [971  has 
been  evaluated  from  Equations  [75]  and  [94]; 
thus 
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The  ratio  K2/Kx  is  given  in  Figure  15  for  a  range  of  diameter  ratios  k  for  viscosity  ratios 
o  *  0,  1,  and  «.  For  the  rigid  case  (o  -  0),  the  ratio  Kt/Kx  obtained  from  the  exact  solution 
Is  also  included. 
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SUMMARY 

The  steady,  slow  motion  of  rigid  and  fluid  spheres  translating  along  the  axis  of  an  in¬ 
finitely  long  cylinder  has  been  examined.  The  solutions  are  given  in  terms  of  the  stream  func¬ 
tion.  An  exact  solution  for  the  motion  of  rigid  spheres  is  obtained  in  terms  of  an  infinite  set 
of  linear  algebraic  equations  for  the  coefficients  in  the  stream  function.  The  drag  of  a  sphere 
in  motion  within  a  moving  liquid  is  shown  to  be  composed  of  two  parts:  namely,  the  drag  due 
to  the  motion  of  the  sphere  in  a  Mill  liquid  inside  the  cylindrical  tube,  and  the  drag  due  to 
the  meti on  of  the  liquid  inside  the  cylindrical  tube  past  a  stationary  sphere: 

Drag  «  6  e  g  /?  (£/  -  K  K2) 

The  drag  experienced  by  rigid  spheres  has  been  determined  for  the  two  special  cases  over  a 
range  of  ratios  of  sphere-to-cylinder  diameter.  A  very  good  approximation  for  the  drag  of  the 
rigid  spheres  is  obtained  by  Utilising  the  first  two  equations  of  tha  infinite  set.  For  fluid 
spheres,  an  approximate  solution  (similar  to  tha  rigid  case)  has  been  obtained. 

The  '^proximate  expressions  for  the  wall  correction  factors  for  the  spheres  in  a  cylin¬ 
drical  tube  take  the  form 

1  -  0.75857  A5 

K  _  _ 1+8/3  o _ 

1  ’  1  +  8/S  a  1  ,  1-8/8  a  .  l-o  . 

1-8.1050  — —  A  +  8.086S  - - A* -1.7068  — -  A* +  0.78608  - A4 

1-fa  l+a  l+a  1+a 


i-l 


8  1+2/8  o 

2~  1  +  2/8  a 

1-2.1060  — -  X  +  2.0665 


X*-  0.20217 


1-a 


1+2/8  a 


1+a 


1  t  1-2/8  a  -  1-a  7 

- —  A*  -  1.7068  — - A*  ♦  0.78608  - - A* 

1+a  1+a  i+a 


Experimental  results  for  tha  rigid  and  fluid  cases  confirm  the  theory.  In  general,  the 
results  show  that  the  wall  affect  for  the  fluid  spheres  is  leas  than  for  corresponding  rigid 
spheres. 


FUTURE  INVESTIGATIONS 

Since  the  experimental  investigation  showed  that  fluid  bodies  deform  iato  spheroids 
(ellipsoids  of  revolution),  work  is  cmreatly  in  progress  to  obtain  a  solution  for  tha  stream 
friction  in  spheroidal  coordinates.  It  is  hoped  that  an  ax  act  solution  for  fluid  bodies  might 
be  obtained  in  this  fashion,  or  it  might  at  least  improve  the  present  approximate  solution  for 
Urge  diameter  ratios. 
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APPENDIX  A 

GEGENBAUER  POLYNOMIALS  c;»U> 

The  Gegenbsuer  polyi.omi.ls  <?;»(*>  *W<*  «PPwr  i  n  t  he  mthturn  tor  the  stre.m  f-etUm 
in  spheric.!  coordin*t*s  c.n  be  ev.lu.ted  by  melting  use  of  their  reLtionship  to  the  Legendre 

polynomials  (Reference  25,  page  77>: 

where  P  U)  is  the  Legendre  polynomial  of  depee  n,  and 
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APPENDIX  a 

EVALUATION  OP  INTEGRALS  CONTAINING  *0(<H>)  AND  Kx(ap) 

The  integrals  containing  K0iop)  were  evaluated  by  snccessive  differentiation  with  re¬ 
spect  to  starting  with  the  known  integral  (Reference  25,  page  8$) 
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-  i  sin2  0(- 10, 895  cos6  0+ 14, 1T5  cos4  0-4725  cos2  0+  225) 

2  r5 

-  —  sin2  0(- 185,185  cos6  0+  218,295  cos4  0-99,225  cos2  0+11,025)  cos  0 

2  r6 

-  -i  sin2  0(2,027,0^1  cos*  0-8,788,780 cos6  0+2,182,950  cos4  0- 
2  r7 

-899,900 cos2  0+ 11,025) 


il 


The  integrals  containing  Kx  Up)  were  evaluated  by  successive  -tsi  nth  re¬ 

spect  to  xy  starting  with  the  known  integral26 


f 


«  f  l 

K.  (op)  sin  ax  da  - - 

1  2  p  r 


Hence 


f,Kta 

f  pKva 

J  o 


cos  ax  da 


2  si  n  ax  da 


I  p  A'j  o3  cos  ax  da 

Jo 

I  pKla,A  sin  ax  da 

Jo 

I  p  K^  a*  cos  ax  da 

Jo 

I  p  k'x  o6  sin  ax  da 

Jo 


f>’ 

J  Q 


^  l  A  W  .  W 

- sin5  9  -  —  CZ*U) 

2  r  f  2 


n  1  ,,  3n  u 

- sin2  0*  3  *  cos  0  *  —  (?»  (() 

2  f2  r5 

ff  1  ,  .  0  12  ir  w 

- -  sin2  0  (-15  cos*  0  +  3)  «  -  - — -  07”  (t) 

2  ,J  f3  4 

- - -  sin2  0  (- 105  cos2  04  45)  cos  0  - - C,  w(t) 

2  ,4  f4  * 

—  —  sin2  9  (#45  con4  9-630  cos2  9+45)  -  — — —  CZ*  (t) 

2  ,3  ,»  * 

-  1  sin2  9  (10, S»5 cos4  9 -#450  cos2  9+  #45)  cos  9 
2  f« 

JT  1 


cos  ax  da  «  —  •—  sin2  0  (- 135+135  cos6  0402,370  cos4  0-36,855  cos2  04045) 

2  r7 

-  sin2  0  (-2,027,025  cos6  04 1,621,620  cos4  0-054,885 cos2  04 

2  f§ 

4  84, 105)  cos  0 


sin  ax  da 


etc. 


Another  group  of  integrals  is  derivable  from 


£ 


1  a?  4  r  1  1 4  cos  0 


A'0(<ip)  sin  ax  da  -  —  In  - -  — -  In 


r  p  2  r  1-cos  0 


ana 


P  *4f  1  P  1/1  1 4 

1  A',  (op)  a  sin  a*  da  - In - — 4 - « - ( —  sin  0  It  — 

Jo  1  ,s  p  rp  f2(l+f)  ,l\i  1- 


cos  0  COS  0 

4 


coe  0  8191  0 


cos  0  \ 
sin  0  / 
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These  integrals,  however,  resuk  in  iyrtaive  velocities  for  p  *  0;  hence  the  integrals 
fmK0  ( ap )  a*  cos  axdc  (where  n  *  odd),  jf*  K0  ( ap )  a*  sin  ax  da  (where  n  *  even), 

rdL  (ap)  a”  cos  a*  da  (where  *  ■*  even),  and  Kt(ap)  a"  sin  ax  da  (where  n  *  odd)  cannot 

o  1  o  1 

be  used  in  ihe  present  problem.  These  integrals  correspond  to  ?m  (t  yterms  in  the  spherical 
coordinate  solution. 


— 


APPENDIX  C 

NUMERICAL  EVALUATION  OF  INTEGRALS  5/  AND  S4m 

The  integrals  S"  and  S4  defined  as 


..  f 

2  J0 

'*  -  f  <<*» 


(«6» 


/,  («4)  K ,  (o4)  +  /2  (a6)  *•„  («4) 

/*(«*)  -  /„  <o4)  /a  (at) 


o"ci: 


I- 


were  evaluated  using  Simpson’s  rule.  Increments  of  0.1  in  the  argument  wer*  used  ?c  :  v3 
range  0  to  3;  for  larger  values  increments  of  0.2  were  used.  v  of  the  Bessel  functions 
were  obtained  from  Reference  27. 

The  S3*  ♦  2S®  had  to  be  evaluated  in  the  combined  form)  namely 


S 


2  (/,  A,  +  l2  K0)  -  1 


/•  2 \L  Kt  ♦  /,  / 

*♦2*®-- 1  - - - 

Jo  '.’-'o' 


ii(ab) 


The  same  increments  as  above  were  used.  The  rt.nge  0  to  0.1  -ov#'  wing  the  small 

value  approximation  of 


i 


o.  1 


[1-2  y  +  2  In  2-2  In  (ah)]  dUft 


where  y  -  Euler's  constant  (0.577215665  . . .  ). 

As  check  on  the  computations  the  integral  S*  +  2S?  was  also  evaluated  in  the  combiner! 
form. 

The  evaluated  integrals  are  tabulated  in  Table  3. 


U 


.jjk  4 


/-s ' 


TABLE  3 


Valuer  of  the  Integrals  S* ,  S",  and  S3  +  2 S  * 


I 


s.+t 


17.87328 


6.55507 


S**l  +  2Sf 

-  4  A  08 66 

-  4.76313 


60.5861 

15.3046 

-  293769 

675.6858 

133.8245 

-  408.0368 

137.7*5$  x  10J 

226.4366  x  10 

-924.7838x  10 

431.9631  x  103 

$09.0)51 x  10J 

-310.1441  <  10s 

1)1.2120  x  10s 

235.514$  x  104 

- 144.1090  x  10s 

112.19$)  x  107 

123.31*0  x  10* 

-  *75.3534  x  10* 

APPENDIX  D 

EVALUATION  OF  THE  DRAG  OF  SPHERES  FROM  THE  STREAM  FUNCTION 


’1 

i 


I 

i 

5 


The  drag  of  the  spheres  is  to  be  evaluated  by  integrating  the  forces  on  the  surface  of 
the  sphere  in  the  ^direction 


D  ■  Dq  +  0 


[D-l] 


where  0  is  the  drag, 

Dq  is  the  tangential  drag,  and 

Dr  is  the  normal  drag. 

From  Figure  16  these  components  are: 


‘■'e  ’ 


J  *r6  sin  $  dS 

Df  «  -  j pfr  cos  0  dS 


where  dS  «  surface  element  «  2  wR2  sin  6  dO. 
Hence 


00-2  »2tr  J 

0, »  -2R*w  | 


pr0a\n2OdO 


prr  sin  Ocoa  OdO 


Figure  16  —  Stress  Components  on  Surface 
of  Sphere 


From  Reference  18,  page  874: 


/i  *«v  *9  5*e\ 

rtQ  « it  I  —  — - —  +  — — - 1  Tsngestial  stress  CD-41 


rr 


-f  +  2p 


l!: 

dr 


Normal  stress 


IO-4] 


Utilizing  the  expressions  for  the  velocity  components  (Equations  [18],  the  tangential  stress 
becomes: 


,4s,  sg  Sv9 
r  BQ  r  dr 


m 


ssU-a^^+j^-DB, 


♦  2n(a* 


...i] 


1 


♦  2(*2-l)  CV*‘  + 


ID41 


IT 


Since  -  (1/2  n~l)  (P_  ,  -  P. ),  the  tangential  drag  becomes 


Dq  «  2/?2  IT 


»— 2  * »' 

*0  2 


«U '  p.> 


|n(n-2)< 


*  '  .  1 

1  R*'3 +  (n2-l)Bm  - + 

Rn±2 


(n2-l)  C  ff*‘l+n(n-2)D  —I 


sin 


From  Reference  25,  page  52 


r 


rr 


P  (cos  6)  sin  6  dO  ~  2  for  m  *  0 


10-71 


0  otherwise ; 


hence 


00-  -^+  C2x\ 


IP-81 


In  order  to  evaluate  the  normal  stress  component  (j>rf),  the  pressure  (p)  must  first  be  obtained 
from  the  equation  of  motion.  From  Reference  18,  page  373,  we  have  for  slow  flow: 

2  8v0 

I?  Tf 


v. 


cot  0 


1  dp  *  y 

—  - -  *  y2  v  -  2  - —  -  2  - -  - 

n  dr  r  .2  r2  0 


[0-9] 


2 


1  dp  2  Va 

- JL  ,  y  v  - - r  —  — 

M  *  r2  sin2  0  r2 


CD-10] 


where  7 


2  d2  2  d  cot  6  d  Id4 


dr 1  *  r  dr 

r2  d  e  *  r 2 

d$2 

d  vr 

1  8  vo 

*hivjp&ia> 

*  —  — . — +  0 

—  + 

dr 

r  dO 

r 

cot# 


2  —  +  — —  ■  0  Equation  of  Continuity  [D-lll 


Hence 


and 


r  \*fr  r  I  \f2  dd  f2  f2  7 


1  /If  2 

— - *  f  r.  ♦ 

H  dr  r 


2  rr\ 

r  \df  +  r  J 


CD-12] 


CD-IS] 


51 


- 


Hence  from  Equation  tt$} 


1 

p  d  t 


27 

2 


»  - 

-1  "S 


2  (2  n  +  1)  C  r*~2  +  2 (2  n  -  3 


D 


*  t-fl 


to-n) 


Integrating 


90 


V  *  f*‘>  ^ 
2 


r  2(2  n  +  n  2(2n-S)  ^  J 

-i  [r  *-l  *  n  "  r* J 


f  f  (0) 


Integrating  the  equation  containing  the  partial  derivative  of  p  with  respect  to  0,  we  obtain 

!"_  C?T*~X  -  2(2n~--  0  -J  +  fir)  [D-16-1 

»-t  ft  - 1  *  »  "  r\J 


2 


hence  /  (r )  *  /  (0)  ~  0  h  —  2a 
Therefore 


pff -p  /^Pm. .  4  -2tn-2)  A.  r"“3  +  2(n+ 1)  S, 
2 


*2 


2  (9  a -8) 

+  i  -  -f  2  in 


/+*  L  "-1  J 

.„],!} 


<7  f*”1  + 

n 


[11*173 


From  Reference  25,  page  52, 

p  sin  2  0rf0  -  0  for  m  >  2  or  m  even 

4 


r 


3 


for  m  <  2  and  at  odd. 


The  normal  drag  becomes: 


y.-lR’ag  (28,  2-  +  2(7,8  + 


[D-181 


Adding  the  two  drag  component,  w.  obUin  for  the  drag  of  a  ephere  moving  in  a  medium: 


0  * -iwuD 


*u2 


The  same  result  has  been  obtained  by  Savic20  for  the  case  A*  ■  Um 


A  -  (7  -  0. 


SI 


•  ‘A 
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APPENOiX  E 

MOTIOH  OF  SPHERES  IN  A  SPlffiRICAL  CONTAINER 


The  motion  of  .  sphere  at  the  instant  it  p.ssee  the  center  of  .  spheric.l  conUiner 
(neglecting  unsteady  effects)*  is  of  interest  as  guide  and  upper  bound  for  the  motion  within  a 
cylindrical  container,  since  the  wall  effects  for  corresponding  spherical  boundaries  exceed 

those  for  cylindrical  ones.** 

The  origin  of  the  coordinates  is  taken  at  the  center  of  the  sphere,  and  the  container 
is  assumed  to  be  moving  at  constant  velocity  V  in  the  negative  ^direction,  as  indicated  m 
Figure  IT.  The  boundary  conditions  for  the  sphere  have  been  given  in  a  previous  section. 
The  boundary  conditions  at  the  container  walls  (f  *  P)  »f«- 

>J  sin2  9 


V 


-  V  cos  0  or  ^  *  ~z  U  P 

£ 


a  U  sin  0 

From  Equations  lF-1],  [39],  and  [^7]  we  obtain- 

A  p*  +  B  — +  0_  P*+2  ♦  P. - r  “ 


[E-13 


0 


*  pn— t 


>a— 3 


1 


n  *  p»+t 


pn 


*  *  2 


[E-2] 

*  0 


Sphere  Passing  Center  of  Spherical  Container 

with  Velocity  V 


Equivalent  Sy atea»:  Sphere  Fined;  Spherical 
Container  Hoeing  with  Velocity  V 


•Tbs  corresponding  two  finaati 
••tee  Reference  29. 


Figure  17  -  Definition  Sketch  for  Sphere  in  Spherical  Container 

•ted  by  Frazer  {Reference  21). 


1  problem  ter  a  rigid  cylinder  bee  bee 


h 


Solving  for  A  ,  B  ,  C  ,  D  ,  £  ,  and  F,  we  obtain 

Jfc  Ji  ~  W  *  *  W 


"V  !J  —.•*> 


1  4  a  - 


5  ,  */«  \  s 

1  4  a  4-  —  X3  -  — {  —  -  a  1  Xs 

♦  g  V  2  / _ 

Ut  +  a)  K+*  X3  -  —(—  -  <r\  Xs  ♦  Cl-rix' 
2 \ 2  )  2  2\2  / 


1  <% 
-  (;£3 

42 


l-(l-a)  X3 


1  £ 
2  ^ 


3 

WJ? 

2 


l  +  <7_l/l  +  „\  X4-X3-i(-~<’)AS><1-'')X< 
2\2  /  2  2  \  2  / 

f|  +  o)  X3  -  —  Xs 

_ _ \«  y  g _ _ _ 

3  / 3  \  5  -  3  /$  \'  *  ' 

1  +  cr.-—  j  —  4.  a  )  X  41  —  X  —  —  /  —  —  a]  X  4-  (1  —  a)  X 

2  \2  J  2  2^2  / 


1  4  —  a  -  (l-a)  X5 
3 


*2~ 


3  3/3  \  5,3/3  \  ,  ,  ,  , 

1  4  a  -  — j  — 4  a  )  X  -4  —  X3 - - 1 - a  j  X5  4-  (l-a)  X1 

2\2  /  2  2\2  j 

„  (i-ix34x*)„ 

&  v  2  2  y 

—  ■  I  ■  !■  !■■■■-- I  ■  -  l«  — — »—  —  ■■!  !!■■■■  —  1  j  l—IIW  — 

2  3/3  \  v  5  ,  3/  3  \  <  ,  v 

1  +  <7'?(2+<’)  +  7  '7(l  ■")  X'+a-’)A 

v  !  (l  -  i  X3  4  l  Xs)  a 

2  Rj  1  +  "  *  7(7 + ")  +  7  x*  ’  ?(!“  0  )  x* +  (1”o)x‘ 


V  1 


A  m  B  m  C  *  D  m  IS  m  F  S  0 

*  A  »  n  w  n 


The  stream  function  in  the  exterior  of  the  sphere  thus  takes  the  form 


*  -  -  sin2  B  U2  r2  4  C2  r4  4  fi2  —  4  02 1  j 


IE.31 


[E-4] 


where  ^2,  fi2,  C2,  and  /?2  are  given  above,  whereas  the  stream  function  in  the  interior 
becomes^ 


(E-51 


whet*  S' j  is  given  in  Equation  [E-S]. 


-*-»*'*  <*»*e 


For  to  infinite  medium  02  m 
for  motion  in  a  spherical  container  becomes 


3  1  +  2/3  a  12 

—  UR  . — „ — —  .  Hence  the  wall  correction  factor 

2  1  +  o 


1- 


-  a 


1+2/8  c 


,  9  1  +  2/8  o  5  1  *  9  l-2/8o  c  1-<t  . 

1 - - - A  + - A3  -  -  — - As  + -  A0 

4  l  +  o  2  l  +  o  4  l+o  l  +  o 


For  a  rigid  sphere  (o  «  0)  the  above  equation  reduces  to; 


H  m 


1  -  A 


9-  5  *  9  -  . 

l>-X  +  -A3-!-A5  +  A< 

4  2  4 


For  a  fluid  sphere  of  vanishing  viscosity  (ft.  «  0  or  o  -  »): 


3  < 
1  +  -  A5 

a 


3  3  ,  , 

1--U  -  A5  -  A6 


2  2 

For  &  fluid  sphere  with  viscosity  equal  to  that  of  the  internal  medium  ( o  *  1): 


[E.7] 


[E-8] 


15  3  ,  3  , 

1 - A  +  —  A3 - A' 

8  4  8 


(E.91  r 


The  wall  correction  factors  for  these  three  cases  are  shown  as  a  function  of  A  in  Table  4  and 
Figure  18. 


TABLE  4  -  Wall  Correction  Factors  for  Spheres  Moving  in  a  Spherical  Container 


Rigid 

Fluid  Spheres 

A 

Spheres 

fh  ’  0 

0.0 

1.000 

UOO 

1.000 

O.i 

1.285 

1.229 

1.176 

0.2 

1.756 

1  *575 

1.421 

0.3 

2.573 

2.126 

1415 

0.4 

4.106 

3.066 

2.469 

0.5 

7.294 

4*831 

3.722 

0.5 

14.948 

8.636 

6.569 

0.7 

37.830 

18.762 

14426 

0.8 

138*224 

58.480 

50.773 

0.9 

1209.778 

431.779 

439.156 

*Thi«  rqvwtion  tgrm  with  that  of  Cun nlnfhwn7  and  U«9nc«pt  that  Caaelughaai  attains  a— A6  tern; 
this  appaars  to  ba  a  print  inf  arror. 


Figure  18  —  Wall  Correct*  «t  Facte  i  for  Spheres  Mox» 


in  a  Safari  cnl  Container 
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